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comme UNDER COMBINED BENDING AND THRUST 


Theodore V. Galambos,1 J.M. ASCE and Robert L. Ketter,2 J.M. ASCE 


__ Interaction curves: relating the axial thrust, applied end | bending moment ; 
al slenderness ratio are developed for the ultimate carrying capacity of pin- — 

3 ended, wide-flange beam- columns. It is assumed that failure is due to — 

cessive bending in the plane of the applied moments which is further con- 

_ sidered to be the plane of the web. The two conditions of loading that are 
investigated are (1) equal end moments applied such that the resulting de- 
formation is one of single curvature, and (2) end moment applied only at one 
SS of the member. ‘The influence of an assumed eyumates residu 


= 


When designing analyzing) a structure the inate plastic theory, 
7 is assumed that the member in question will deliver the fully | plastic moment — 
value, Mp, noted in the calculations. This, however, will not necessarily be ae 
the case if the member is subjected t to an axial thrust in addition tobending 
_ moments. (1) To attain the desired moment value, itis necessary to supply : a Re 


"without modification to rolled wide- columns of A-7 structural steel. 


_ Similar procedures could be be used for other her materials }and loading conditions. 


te: Discussion open until September 1, 1959. To extend the closing date one month, — 
ie written request must be filed with the Executive Secretary, ASCE. Paper 1990 is 

_ part of the copyrighted Journal of the Engineering Mechanics Division, —— 
of the American of Civil 85, No. EM 2, April,1959. 
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Bt = _ for slenderness ratios up to and including L/; = 120. For ease of future © 
— 
< el pure bending; i.e., one that will develop the required end moment inthe = a * 
Prot. and Head, Dept. of Civ. Eng N. Y.; former= 
4 ly of Fritz Eng. Lab., Lehigh Uni 


However, to: afford ; a means of with predicted 
= and to facilitate solution of problems where the material has a = a &s 
- yield point other than 33 ,000 psi, a method of modifying the solution is me ; 
The ‘problem that will be considered in this paper is the determination of 
the maximum amount of end bending moment that a member can sustain when J : 
it is also subjected to a given axial thrust; the material presented herein con- - . 


_ stitutes an extension of certain of the ideas advanced in an earlier paper.(2)_ ‘7 
Two loading cases will be investigated: 


* axial thrust plus equal end m moments applied at both ends ies member = 
such that it deforms in single curvature, and 


2. axial thrust plus moment applied only at one end of the member. sat 
d 


a 


Fie, ‘These conditions are shown diagrammatically as loading conditions “ c” an 


‘ 


Fig | CONDITIONS OF Lone 


ideally -plastic; i.e., there is initially a linear range wherein 
= is followed by a constant stress level a= Oy for strains greater r than 


a €y.* This type of behavior is typical of mild structural (ASTM A7) steel if a 

i _ strain- ‘hardening is neglected. There is, however, assumed to be a symmet- 
oer vical residual stress pattern present in the member prior to the rh rong 
5 a of any external loads. The presumed pattern (Fig. 2) is consistent with 

measured residual stresses in wide-flange column type sections resulting 
_ from cooling of the section during and after rolling. (2 
As shown in Ref. 2, if the material is homogeneous and isotropic and w... 
bending strains are assumed to be proportional to the distance from the neu- — 

oe tral axis, then the thrust- moment-curvature relationship for the 8WF31 sec Bi 

tion will be that given in Fig. 3. Ih this figure two conditions are illustrated. “ 

one lines are for the cases where residual stresses are neglected. 


1s given in Section vill. 
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Fig.2 ASSUMED COOLING 


Sten sla “STRESS” PATTERN 
Toe solutions which include the influence of the residual stress heehee 
shown in Fig. 2 are given by the dashed lines in Fig.3. 
Be. Since the basic approach that will be used in solving the problem con- 
a _ sidered in this paper is one of numerical integration, and since this integra- 


will proceed from a knowledge of the curvature values of Fig. 3, 
in- 


‘tt should be noted, however, that this Section has one of the more severe 


thrust- -moment-curvature relationships of the column sections rolled because > 


of its low shapefactor (1.11 as compared to 1.14 for most sections.) Using 
_ the interaction curves for other shapes should therefore result in a some- 
_ what conservative or at least equal prediction of strength for the member in 
= For ease of presentation and generalization, load and section property 


30 ,000 psi. Since specifications require a minimum yield of foy= 


= 33 ,000 psi for A7 steels, this — was used in the calculations as the bale . 


The approach that will be used in the solution of the problem in /question y 
"will be one of numerical integration. (4) | The calculation will proceed from % 
assumed deflection configuration and will take into account the non- linearity a 
between moment and curvature as strains exceed the initial yield strain. ce 
«Since deflections must be assumed, it is desirable to know the equation of 
2 the column centerline at initiation of yielding for each of the conditions of aa 


loading. These can terme of the from a consideration ¢ of the on 


a are as follows: 
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THRUST- CURVATURE 


RELATIONSHIP, INCLUDING 


INFLUENCE OF RESIDUAL 


4. 
be a) Moments applied at both ends of the member, (condition “c” in Ref. 2) 
ition “d” in Ref. 
For the assumed values of E = 30, ,000, ,000 pst and Oy = = 33 000 2 
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. = Equations (1) and (2) it can be seen that for the conditions of con- 
- stant axial thrust and elastic behavior | there is a linear relationship between — " 
+ the applied end moment, Mo, and the resulting deformation. The maximum — 

‘@ _ Value o of Mo for which this situation holds is referred to as the initial yield | 
g value and the solution to this problem has been presented in Ref. 6 and else- 
where. For greater values of applied end moment yielding will occur at the 
most highly strained sections along the member. . In these regions the mem-_ 
ber becomes relatively weaker to further increases in loading. This can be a 
os from the moment-curvature diagrams of Fig. 3. The load-deformation 
‘relationship of the member as a whole will also indicate this decrease . a : 


_ the fact that the total deformation is the integrated effect of all of the —_— 


- values along the length ofthe member, el 


“a To be able to determine the maximum carrying capacity of a given mem 
ber, i it is essential that the load-deformation of that particular 


the member and successively correct these assumptions based on 


=a and axial thrust ratio, then, the definition of the load- -deformation | 
relationship above the elastic limit, and thereby the definition of the critical © 
loading, may require the consideration of four or five end moment values 
which in turn may require three or four numerical integrations each, tw 
In addition, for a given slenderness ratio, itis necessary todetermine 
s: critical value of the end moment for various values of the axial thrust. 
‘This would make it possible to define the relationship between axial thrust a ; 
a and end moment for this one particular | slenderness value; i.e., to define the a. 
interaction curve for this given slenderness ratio. In general 0.2 Po/Py in- 
tervals were used in the computations on which the interaction curves of this — 
report are based. For a better definition of the ‘relationship | at higher values _ 
of thrust, however, a closer spacing of values of Po/Py had to be used. vita 7 


Slenderness ratios ranging from 0 to 120 were considered in intervals of 20. _- 


In outline form, then, the steps that were used in determining each of the — 


interaction curves presented in this report are as fpllewa: 


GIVEN: loading condition, slenderness ratio axial thrust 
value, for the 8WF31 Section used as a standard. 


1. Assume an end moment, Mo, greater tl than the initial yield value; 


"assumption of deflection would be satisfactory. 
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2 Computing the moment values a at eight Spaced s stations 
along the length of the member (Mx = Mg + Poy) and numerically a * 
integrate curvature values obtained from Fig. 3 (an enlarged ver- 
_ sion of this figure was used). , Figure 5 shows cme cycle Fe 


4. Correct the assumed deflections based on the values obtained — 
from this numerical integration and repeat step (3); 


step (4) until the desired accuracy is obtained 0.001 
was used in this report); = 
critical 


no 
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6. the end rotation the final deflection values of 


if it is 3 assumed that the deflection curve of the member within the three ce 
expressed in terms of the known deflection as" 


deflection at first station away from the end of 
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RESIDUAL STRESS) 
Assume values of the | moment, Mo; and repeat 
8. Plot the various values of My versus 66 from step (7) and deter- ; 
mine the maximum value of Mo the curve. . Gee 


Se deflection at the second station away from the applied a 
Az grid spacing (assumed to be L/g for the cases considered) 
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Fig 7 MAXIMUM CARRYING C 
INFLUENCE OF RESIDUAL STRESS) 


This gives one particular point on one suai interaction curve. As w was 
pointed o out above, ‘it is necessary toc determine r= points to be able 

define the desired range of the interaction curves. A 


the ne (Mo/My)critical values obtained from the numerically deter- 
Za Mo versus 85 curves by the shape- -factor (i.e. ey ), the interaction bl a, 
curves Po/Py versus Mo/Mp versus shown in Figs. 6a and a7 were ob- 
tained. Fig. 6 is for the case of moments applied at both ends ¢ of the ‘member 
D (condition c”) and Fig. 7 is for the case of momert applied at one end (con- ce ' 
Only the interaction curves incorporating the influence of residu- 
al stress have been included in this report. However, interaction curves \ ld 
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ial 
; , a tion curves are shown in Ref. 7. Also given therein is a more detailed expla- Pi -— 
nation of the derivation of the curves shown in Figs. 6 and 7. To give an ce J 
indication of the influence of residual stress on the carrying capacity of 4 
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RESS 
members of the type herein, Fig. gives 
— make the curves more useful when eccentricity ratios (ec/r“) are given 
ars 
nana of end moments, values of ec/r2 are also shown on the interaction th 
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BENDING AND 


To avoid from om the of 6 1, approxi 

ae interaction equations | were developed by fitting the curves into c cubic | and — 
quadratic equations. All of the limitations of the original curves are there- 
fore present in these approximations. In general, the range of application an 

_ was chosen as 0 < L/y ¢ 120 and0 ¢ Po/Py 0.6. It was anos oa 


er. 
these covered the 1 major range of practical applications. 

1. Pin-ended column subjected to axial thrust pilus two equal end-moments — 


lied such the deformation is that single 


Assuming an the form 


‘where K and are to be the slenderness the 
coefficients of the axial load terms of the “ery vers found to be 


-0.a20 +H) — | 


29, 00 


©. 770 — 


= agreement between these expressions ao the « comparable relationship: 


tion between eskess obtained from the approximate th (shown a as 


points) and the interaction curves of Fig. 6 is given in Fig. 10. oO. 


i. -ended column subjected to axial thrust p lus an end-moment lied 
only at one end of the member (condition *d”): 


yes an equation of the form ihe; 


o. 
‘Py 
_— oe in case 1) B and G are assumed to ve: functions o1 oniy of the — fi 
ness ratio, the coefficients are found to be 


a 
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ACTION ‘CURVE OF Fig6 
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ATION 
SOLUTION OBTAINED BY NUMERI- 
CAL INTEGRATION (Fig. 


BETWEEN EXACT 
ACTION CURVES ( LOADING 


It be noted that when Equation (7) predicts a of 
1.00 (that is, for small values: of Po/Py ), Mo/Mp = = 1. -00 should be 


padre determined numerically (Fig. 7) is shown in Fig. 11. —— 
Table 1 is a tabulation of the interaction equation « B, G, Jan 
L/y values from 0 to 120 in 


an 
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SOLUTION 


COMPARISON BETWEEN "EXACT" 


"APPROXIMATE" INTER- 
ACTION CURVES (LOADING CON 


DITION"d" 
= maximum axial thrust that the member will sustain ven 
gubjected to pure axial thrust; 
 Mt=m maximum end | moment that the member will sustain 


* Euler buckling load for the axially loaded member. 


it assumed in the derivation of the interaction curves presented 
earlier in this report that the member did not fail by lateral instability, M' a 
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* Axially Loaded ed Columns 


645 
a shown in Figs. 6 and 7 where M,/My = = 0, a comparison has not been shown. © 
It is gratifying to note the close Correspondence between Equation (11) of — 


_ this paper and Equation (20) of Ref. 17 by Bijlaard, Fisher and Winter, since 
— the latter expression was determined by an entirely different —_— - 


Asan experimental check of the theoretical predictions of 
- existing test data are compared with the interaction curves of Figs. 6 and 7. i 
|The tests of the following experimental are used for comparison: 
Lehigh University, 1940 (Ref.10) 
3. University of Liege, 1956 (Refs. 11, 12) 


4 University of Wisconsin, 1920’s (Ref. 13) 


_ ‘Lehigh University, , current 


_ Graphs comparing the analytical geen with experimental r results a are 
shown in Figs. 12, 13, 14, 15 and 16. 


i ne ‘In Fig. 9, the e ssions for K' and J' (determined from Equations . 
ME §=—=—=—séing for (P* /Py) the end points for the case where ;,° = 0 of Fig. 6 « a 
on ip between axial thrust and 4 4 
slenderness ratio that occurs below the Euler curve (that is, 0 L/, £112), | a 
— = 
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| COLUMN TEST _RESULTS(!3) 


Equation (12) gives the adjusted slenderness ratio be used 


COMPARISON | OF WISCONSIN ft 


33 ,000 psi with the predictions of pre herein. In the following 
r comparisons the slenderness ratios have been adjusted according to Equa- — 
~ tion (12). Since the axial thrust and end moment parameters were given in a 
Be r ‘non-dimensional form involving the yield point stress, it is not necessary to 
_ For a majority of the columns that have been tested and are listed herein, 
a the members were subjected to eccentrically applied thrusts. In ner 
_ comparing these test results with the strength predictions of this report __ 
_ Slenderness values have been shown as the and values 
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6 using the the values shown across the top along the right 
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re and since e the values of ec/12 given in Figs. 6 and TV were obtained | by using 
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which fully prevented lateral - torsional buckling and therefore - 

_ formed to the assumptions of this report. Two “Z” sections were welded 
together in the form of a “hat” by intermittent welds. Bending was forced iy 
(by the use of knife- -edges) about the minor axis of the total cross-section. 

_ - Figure (12) shows the comparison between the test results of Mason, | 
Fisher and Winter and the theoretical predictions of Fig. 6. In general the 
correlation is quite good. The experimental r results fall slightly above the 7 
_ predicted curves as would | be expected since the shape factors of the sec- 

— tested (f = 1. 18, 1.25 and 1. -17) were greater than those of the 8WF31 
"section. Also, the residual stress distributions of the sections tested were 
‘not as severe as those that were assumed. 

La _ Table (2) gives a tabulation of the data from which the test points of 
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— mea , 
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— ‘ted Beye 57-6. = 


a BENDING AND THRUST 
- A total of 93 tests were — ‘out by Johnston and Cheney in this Ot 
89 were made on 315. 7 sections and 6 on 6WF20 : A summary of 
the test dataisgivenin Table (3), 
Columns were tested by both conc -entric and eccentric application of the 
axial load; however, the column tests under pure axial thrust cannot be com- 
- pared with the predicted interaction curves. The end-condition of these test 5 
. specimens were such that they fail by buckling about the weak axis. 2 
ae _ In general, the tests were performed on columns which were essentially 
pin- ended with respect to bending in the strong direction and fixed-ended — 
the weak direction. This was accomplished by the use of knife-edges placed ; 
Perpendicular the web through which the load was applied. The loading wal 
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the yield stress of the | ‘material tested. The comparisons between eta: . ] 


_ Johnston and Cheney report that the “columns loaded eccentrically to pro- - ; 


duce bending in the strong direction usually failed by plastic lateral torsion- 
_ al buckling” (a condition specifically excluded in this paper). It is interest- — 
ing to note, however, that except for the tests which fall close to the case — 
where failure would have been due to Euler buckling in the weak direction, 
the correlation achieved with the — theory which neglects breton, 


of University 0 y of Liege ge (Ref. 4, 12) ‘wld 

Massonnet reports the results of 95 column te tests were 


- out on sections of DIE 10, DIE 20 and PN 22 profiles. Of these, the DIE pro- a 


_ being considered in this report. | Therefore, only the DIE profile tests will — 
_ be used for comparison. _ Furthermore, only those tests which correspond to 


files are geometrically similar to the American wide- -flange shape, the shape ‘ 


; _ The end conditions for Massonnet’s test columns were essentially pin-- 


the condition and “d” loading are listed. q 


ended in both directions since the end-fixtures consisted of almost friction- 
hydraulically seated steel hemispheres. For such end- conditions, the 


Be least possible restraint is provided against lateral torsional buckling. al 4 


_ Table 4 summarizes the applicable test data for the DIE profile tests. = 


7 Figure (14) gives the comparison between the tests on members subjected ih 
a a condition loading and the theoretical predictions shown by the dot- 


curves. As before, the slenderness ratio is adjusted for differences in yield 
oe level. In all cases the final failure was by lateral torsional buckling; 
in spite of this, most of the test points agree eather well with the theoretical 
relationship that neglects this type of failure. It is expected that further — 
Ee work taking this mode of failure into account will provide | a better 
understanding of the aise and will es in a better correlation in the _ 

comparison has been condition “c” tests of Ref. 12 

- (equal and opposite end moments). Due to the condition of loading and end 
restraints it would be expected that lateral-torsional instability would a 
prior to the theoretical load predicted in this paper and this was found to be- 
the case. A solution to the problem of lateral- -torsional buckling for this 
condition of loading which also takes into account the influence of residual — os 

_ stress has just been completed. | 

report latter ‘soon be available, 


The members tested in this investigation were 8H32 shapes, similar t to the 
at currently available 8WF31 section. The end- -conditions were essentially pin- 


7 a ended against strong axis bending and fixed in the weak direction. Of the five 
> ae tests carried-out, the two which had an adjusted slenderness value greater | 
than 50 failed by lat lateral- ~torsional nal buckling a ata load less 


conditions and support arrangements for the tests correspond to the condi- 
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Lehigh University (Ref. 2, 6, 16) 


‘Table (6) summarizes the results ( of the tests in this series that are 


ric 


most of these have not ‘been shown on the pure axial load 
however, Fig. 16 shows the correlation with predicted strength. An addition- om 
test by (4WF13 » L/r=130) ha has been to extend the range 
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Vv. "DISCUSSION AND SUMMARY 


oe For the conditions of end restraint and the loading conditions he Fig. 1, 
g solutions to 9 the problem o of the determination of the maximum carrying ant 4 
capacity of wide flange shapes loading in the plane of the web have been pre- a 
4 _ Sented. . These solutions assume that the member in question will fail by ex- 
cessive bending in the | plane of the applied moment. Failure due to. lateral- @ 
5 torsional or local buckling has not been considered. The resulting interaction i 
curves | (Figs. 6 and do, however, include the influence of a i 
Approximate interaction equations, which cover the range most prey en- 
countered in practice, were developed 1 to eliminate the need for interpola- bg 


Currently available test results were against | the strength | pre- 


were the only o ones that directly fulfill ‘the assumptions ‘of this ‘report and the 
_ correlation was shown to be very good (Fig. 12). For the cases where the E 
8 members tested were pin-ended in the strong direction and fixed in the weak, 
r the curves give reasonably reliable results provided Euler buckling inthe __ 
weak direction was not imminent (Figs. 13, 15 and 16). Where the members 
_ were pin-ended in both directions, lateral-torsional buckling was a major 2 | 
faa factor in determining strength (Fig. 14). The test results corresponding to a 
ae this situation (Massonnet) seem to indicate that for a condition “d” type of _ 
a loading the overall behavior can still be approximated by the curves of Fig. 
(7) but with less accuracy than in the aforementioned cases (Fig. 14). Mem- 
_ bers loaded in a condition “c” manner, however, carry ——- less load at 
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that were > held constant are e underlined. 
Oy = determined from tension cot coupon tests. ple 


# Adjusted oy (to take : the influence of rat) 
peal was obtained by pro-rating the tension coupon value in pe, 
= same ratios as those given in Ref.15. (Note: values 


-Purther work is currently underway to include the influence of 
_ torsional instability into the strength calculations } and preliminary results ae 
of this study indicate that good correlation = be achieved when this } type of . 
failure is considered. 


: 


4 Frames and Their Components” currently being carried- -out at Fritz 
: ‘7 neering Laboratory, Lehigh University under the general direction of Lynn 8 s. 2 
‘Beedle. The investigation is sponsored jointly by the Welding Research = 
4 Council and the Department of the Navy with funds furnished by the oping 
Iron and Steel Institute, American Institute of Steel Construction, Caine © of “i 
Naval Research, “Bureau of Ships and Bureau of Yards and Docks. 
bs William J. Eney is Director of Fritz Engineering Laboratory and Head a. 
the Department of C Civil Engineering. 


— 
— 
— 
> 22 | 0.82 | 0.23 | 0.887 | 
— 0.12 | 0.76 | 9.12] 1 | 
rs 
= 


April, 1950000 
REFERENC 
DA 


Beedle, 
_STRUC "TURAL STEEL- LECTURE NOTES,” AISC and Lehigh University, 


2, Ketter, R. L., Kaminsky, E. L., and Beedle, L. Be: “PLASTIC rn 
_ TION OF WIDE- FLANGE BEAM Transactions of the 


4 ” STRENGTH OF STRUC TURAL MEMBERS,” Proceedings of the Tth 


Session of the Research Council, 1957. 


Newmark, N.M., “NUMERICAL ‘PROCEDURES FOR: COMPUTING DE- 
FLECTIONS, MOMENTS AND BUCKLING LOADS,” of the 


_ ASCE, Vol. 108, 1943, p. 


Timoshenko, S., “THEORY OF ELASTIC STABILITY,” McGraw- 


and Johnston, B.G., “COLUMN STRENGTH 


q UNDER COMBINED BENDING AND THRUST,” The Welding Journal © 
31(12), Research Supplement, 607-s to 662- 1952. 
(12), 


Galambos, V., and Ketter, R. In, «FURTHER: STUDIES OF COLUMNS) 
UNDER COMBINED BENDING AND Fritz Engineering Labor- 


atory F Report N No. 205A. 19, June 19! 1957, 


A Report of Committee E Research Council “SOME 


RECOMMENDATIONS RELATING TO DESIGN SPECIFICATIONS FOR 
STEEL BEAMS AND MEMBERS SUBJECTED TO COMPRES AND 
BENDING,” May 1954. 4 
8. “Focm R. E., Fisher, G. P., and Winter, G., “TESTS AND ANALYSIS OF al 


| SECTION, Progress Report No. 2, American Institute of Steel 
C., and Campus, F., “CORRESPONDENCE ON 
‘STAUCHION PROBLEM IN FRAME STRUCTURES DESIGNED ACCORD- - 
= TO ULTIMATE CARRYING CAPACITY’ by Horne, M. ,»” Proc. 


Inst. of Civil Engrs. | par Il, Vol. 5, pe p. 558-571, August 1956. 
nm Massonnet, C., and Campus, F., “RECHERCHES SUR LE FLAMBENTENT 

on DE COLONNES EN ACIER A 37, A PROFIL EN DOUBLE TE 


OBLIQUEMENT,” I.R. LA. Bulletin No. 17, April 1956. 


1. Gozum, A. T., and Huber, A. W., “MATERIAL PROPERTIES, — 7 


STRESSES AND COLUMN STRENGTH,” Fritz Report N 
220A.14, May 1955 (Revised November 


iy 
7 = 

— 
= 
| 
4 
OL Eng., Cornell Univ., ithaca, New York, April 

= 
| 

— 
COLUMNS” Paper No. 1789, ASCE Transactions, Vol. 95, 1931. 
Huber, A. W., “THE INFLUENCE OF RESIDUAL STRESS ON THE 

4 
— 


BENDING AND AND ) THRUST 


16. S., Ready, J. A., and Johnston, B. “TESTS OF COLUMNS 


oe ie UNDER COMBINED THRUST AND MOMENT,” Proceeding, g, Society for 
Vili, No. 1 1, 109 (1950). 
“ECCENTRICALLY LOADED, | 


Bijlaard, Pp P., Fisher, G. D., Winter, 
END-RESTRAIN ED CO COLUMNS” Transactions of the ASCE E, Vol. 120, 


18. Shanley, F. R., “APPLIED COLUMN ‘THEORY? Transactions of the Decs 


Area of cross-section (in m2) 


N on-dimensional constants 


-Young’s Modulus ¢ of Elasticity 000, p00, psi 


Length ¢ of member 
Bending of Moment (inch- -kips) 


moment at the end of the 

‘Fully plastic moment value under pure moment 

ar Initial yield moment value under pure moment — 

Axial thrust kips Wo 
Axial thrust at load capacity for beam- 


‘= ection modulus about the strong a axis 


a, Plastic modulus about the strong axis (in 
5) 


of gyration about the strong rong axis 


Thickness of web 


1 | 
4 
a 
— 
M)=Z0y — 
— 
— 
> 


Slenderness rat ratio. 


Non- -dimensional 
Constant defini roperties of material * 


a Ve 
Curvature corresponding to initial yield under pure 
Length of spaced segments of total 


“Yield stress to be 33 k ksi for AT steel) 


20 org Distance along the axis of a member, as shown on 
| 
— 
— 
en 
y Ed | 
— 


NGINE 


_ SEEPAGE LOSSES FROM IRRIG 


. Hammad! 
4 


This paper deals with the two- | of seepage flow 
under gravity from a canal into a. semi-pervious clay layer of finite thickness — 
7 underlain by a freely permeable layer of sand and gravel in which the piez- e 

- ometric head is very near the canal water r level. Two steps of conformal ; 

_ mapping are used and an approximation to the canal profile is adopted. 

this approximation the specified canal width and 


many agricultural zones, especially | in river the stratified 
is formed of a | top clay layer of finite thickness and medium or low permea- ee, 
- bility, underlain by highly permeable layers of sand and gravel. The piezo- — 
metric head in the lower freely permeable layers is controlled by the power- 


soll 


"the top clay layer are usually run at higher water level than the piezometric _ 
level in the lower sand and gravel. Consequently a continuous steady Seepage — 


flow takes place . from the canal through the clay layer to the sand and gravel tp ; 


met with in places with high ground water-table. The same problem, lh 

‘ with very deep lying water- table in the lower sand and gravel layer has been ti 
solved in 1934 by V. Vv. Wedernikow. lo 


_ Asm may be found in texts « on flow through pc porous media, (1) the governing | 
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i s be satisfied everywhere in the field 0 ? b 
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where k is the of the porous mé¢ medium, is the pressure at 
= P is the fluid density and g is the acceleration due to gravity. ee 


a to the above assumptions, the b boundary conditions are: 
‘The line of separation the semi-pervious clay strip and the 


lower freely permeable layer of sand and gravel is an equipotential line. — 


2. The « canal is also | an an equipotential line. 


3. The pressure vanishes ¢ on the free surface which is 

4, The rise of water i capillarity from the free water eutees is 


Method of Solution 
ahaa first order solution of the problem of seepage from a a shallow canal — 
(depth d < semi-breadth B) may be ontained by considering the seepage under 
through an n aperture (-b<x + +b) in the solid boundary y = 0, into 
_ the semi- pervious strip bounded by | y = 0 and y=Din the Z- plane as shown aos 
in figure. The line segment vy = 0, -b< x< +b) is an equipotential, while = 
= =0,b< x < co) ) and (vy =0, <x <-b) are stream-lines. The line D q 
: It has been already assumed that the free water surface in the ern 
layer is high and very near the canal water r level. 
= approximate solution, the line segments (y = °, ,b< x< oo) and ( y = 
=00 <x < - b) have been taken as the top stream -lines instead of the real 
- free surface which is very near to them, as our original assumption. hella 
_ _ The complex- potential w (= ? + ip) of the flow eit be obtained by two 


“each step, being the velocity - potential, the stream-function and i i= 


functions: 
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the points having the same numbers pss mah ia correspond to one anotl 
and A isa constant. Integrating relation (3) we get 


| 


— 


7 

gives the ‘required transformation. 

- 2 Following | the same steps in mapping aia the rectangle 
bounded by y= — K, ¢ = a. in the w-plane onto the lower half of the 
(see figure), we finally get in terms of Jacobi’s elliptic 


_ where k (= tank b/2D) is the modulus, and K, K" are the complete elliptic © 

integrals of the first kind corresponding to the modulus k ; k and the eal 
mentary modulus (= 4/1 - k2) respectively, 

a From (6) and (1) the w-z relation of the flowis 


For a given c: canal of breadth 2B - depth d, an the 
aperture and passing , through the points Z=t B, Z = id, is a fair approxima- 
tion to the canal profiles met with in practice, 

"where - -c is along the canal profile. Fora specified breadth : 
and depth d of the canal, equations (9) and (10) give the two unknowns c and ‘, 
. where | k= = tanh (xb/ 2D). In other words the breadth and depth of the canal | 
define the breadth of the aperture 2b and the constant poten- 
tial value to be assessed to the canal profile. Let B' and d' 
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n(e#) 
(12) gi give on solution 


~ Length of the canal normal to the yee of motion is accordingly given = a 


paal 


where h is ‘the hea head loss between the and the | line y= D or the 
‘difference between the water level in the canal and the free water surface in = “fe 
the soil at great distances = 
_ ‘The seepage flux given by (14) is slightly exaggerated by a small amount . 
_ corresponding to the flow in the narrow strip between the real free surface 
and the approximated free surface (y = 0) which is assumed very neartoit. _ 
| may be worthwhile mentioning that in the case of a drain in which the i 
water level is lower by an amount h than the ‘piezometric. head in the = 
and gravel aquifer, the flux will be reversed and the drain will be gaining — 
_ water at the rate given by (14) instead of | loosing. 7 eth ale ee 
From relations (2) and (8), the equation of the free water surface may be 


For deep canals @ a similar approximate solution may be obtained 


_ciding with the y-axis and the remaining boundary conditions are as 


For small of semi-circular a point | source solution 
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4 layer of thickness D D = 50 ft. and sn Sebbincibnite k = 
the piezometric in the sand and gravel a onal 
= tan (na/20) 2 x 5 
= tanh (B/2D) = ‘tanh )= 
Substituting in (13) ) we get 


‘a the modulus of the elliptic: integral k k = - 0.717, and | the potential ; along the a 


= 


canal profile -c= -.220. gall ty ‘wort 


2 


“ From tables of elliptic fi functions, s,(4) tt the ‘complete elliptic integral of modu~ 
lu 


us k = 0.717 is K = 1. 860. . For the complementary modulus k' = v (1 = k2) = 
0. 695, the complete elliptic integral K'=1.835. 
hd Hence the e seepage | loss } per | unit nit Length ¢ of canal | as given by (14) )is- 


011 gal./min./ft length of 
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7 TORQUE- -LOADED CONTINUOUS BEAMS ( OF PROFILE a 
D. H. and Brahtz, M. ASCE 


SYNOPSIS. 
of a beam of thin profile. section is continuous over several supports. and 
subjected to transverse loads not on the shear-center axis, it will undergo — 7 
——- as well as bending. Beginning with the general theory of non-uniform - 
torsion for such bars, expressions for the angle of twist ¢ along the bar are 
_ written for two consecutive spans. From these expressions, a single recur-_ 
rence type formula is derived which involves, as unknowns, only the three — 
quantities at three consecutive supports. ‘This recurrence 
formula is somantat analogous to the three-moment equation and greatly 
_ simplifies the analysis of the torsional behavior of a continuous beam sub- 
jected to eccentrically applied transverse loads. Its application to the case a 
of a two-span beam is illustrated by an example and an experimental study ae a 
of this two-span case is described. The results of this Pees: are in 


- ae Fig. 1, the case of a continuous beam of thin-walled | open 


Se section subjected to the action of various transverse loads. gl 7 


“paper is to develop a general method of finding the stresses a and deformations 
by the torques. can be done by using the general theory of 


non-uniform torsion of bars of thin walled | open section as developed by 
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a torque. The loads on the Shear-center axis produce only bending and 
e corresponding stresses and deflections can be found without difficulty 
. - 


poilew 


tr 
= 


‘aan (2,3) Then the differential equation — the relation —— 
T and angle of twist ¢ along the bar is as 


_ wherein mii primes denote differentiation with respect to z. The torque T is 
considered positive as shown in Fig. 1 and C, Cj, respectively, are the tor-_ 
; = rigidity and warping rigidity constants for the cross-section (1). —- - 
Consider now any span, Fig. 2, a concentrated Pe 


support. rt. For su such case, the solution of Eq. (1) ) becomes 


where = -C/Cy and 


a 
q 


When these expressions are for A, in of 


“which ses the angle of twist: at any cross- z in of 


Fora at zz {, the general for can be obtained by 
posing on the eolation a a term to account for the effect of the 


+ £0 sinh kz + ke ~ 1) - 2M einn kz - - ke) 


aii (2) and (3) represent the complete solution for the case of a > al 
span subjected to a concentrated torque Pe at some intermediate cross- | 

a all If, instead of a concentrated load P, a single span is subjected to uniformly _ 

, distributed load of intensity q applied with eccentricity e, Fig. 2b, the torque & 


Ps To - - qez, where the at the left ¢ end of the taken as 


Timoshenko,(1) Wagner,(4) and others. Such treatment willbe much simpli- 
a 
id ie ofthe left-hand 
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a constants of integration. For z= 0, this 
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inh ke + + kz - 1) - 


(cosh kz - -1) 7 


_ Consider now any span {,, of the continuous beam in Fig. 1 where it is as-_ a 


_ _ sumed that all supports are on the same level and offer no more a 
fe 7 _ than to to prevent x- and y- displacements as well as rotation about the Z- axis 
of the ‘supported cross-section. n. Then at the ends of this span, on = = = 
Substitution ~ these end conditions together with z= hh into Eq. (3) and its 
derivatives and gives, , for the case of. of a a single 


(Pe)n on the span, 
wi rrr 


on kf, + 6" cosh 


Elimination of ian between the first and third of these e equations yields <a 


just to the right of the support 
of g'' and ¢" at the two ends “of the span {_ and the applied torque 

load (Pe), ne There is one such equation for each span of the continuous beam. 

~ 

torques Ty just to the right of each ‘support, can be computed 


from Eq. (6a) and the problem is solved 


To find a relationship between three consecutive values 


three supports, Eq. (3) can also be to the ‘span This will Lave 


q 
set of three e ‘equations to Eqs. (5) ‘above, ‘except that the 
_ Script n becomes n+1 while the subscript n+1 ‘becomes n+2. _ These Six equa- i 
tions for two spans contain eight unknowns, namely, 


tedious, it is 


Py ol a a the origin of z. With this value of T, integration of Eq. (1) gives eT ae a 7 
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where the functions E,”,W, depend only on the physical di dimensions ofthe 
ane the of the functions are defined as follows: 


¥ 


“> 


i Equation pe is a recurrence a — which may be written once for 
each pair of consecutive spans throughout the continuous beam. Thusits 
application t to a beam subjected to various applied torque loads is somewhat _ 
_ analogous to that of the well known Three Moment Equation used in the analy- 
sis of continuous beams subjected to bending loads. In the case ofabeam ~~ 
“on N N+1 supports and having N spans, Eq. (7a) can be written N-1 times. 3. This 
‘gives a of N-1 linear. equations involving the N+1 unknowns 


on 


if the cross- -sections at the two extreme ends of the ’ 


are free to warp, N 470 and Eqs. (7a) are sufficient to deter- 
eee oie As soon as these values are known, the internal 
_ torques Ty Ter Ty _ can be found from Eq. (6a) and the p 

ll the neal left end of the beam is built-in so that warping of li aa 


cross-section is prevented, the value of oy is not known but instead, as a 


‘= 0 vill b b Ss ti i a. 
n w e sa Ea. ‘be for ‘these 


anil stituted for both {pn and {n41. The same procedure may be used if the ex- 
_ treme right end of the beam is ee in. ic hea 
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flects the span {1 and its torque load (Pe); as shown in Fig: 3. Then the re- 4 
_ flected span {1 and the actual span represent two consecutive spans for __ 
— 
con yULLVe Spans yie ng la 10nal equation to matc ea i0nal } 
sits 
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1950 
o The previous discussion has treated ity the case of a continuous beam ae 
. subjected to eccentrically applied concentrated forces. The case of uniform- 
ly distributed load of intensity q applied with eccentricity e from the » Shear~ 
center axis, Fig. 2b, can be { treated in the same way, using the general solu-— ‘2 
tion represented by Eq. (4). _ this case, the recurrence formula, analogous 


= 
. where bag) and (kf) are the same functions already introduced in Eqs. | 
and th function ¢(kf) is = 
- e new function 1 S(kf) i 


| 


case of a continuous beam under the action of uniformly distributed load. ui —% 


_ Tf there are several concentrated loads in any span, or both concentrated ¥$ ‘ 


appropriate right-hand members of Eqs. (7a) and (7b) and likewise ” | 

terms in Eqs. (6a) & (6b), 

+a To facilitate the use of the recurrence formula (7a) or (7b), the es “= 
é,n, ¢, and y, defined by Eqs. (8a) and (8b), ‘are recorded in Table I for a _ 


- loads and uniformly distributed load, it is only necessary to superimpose — 


Experimental 


ponding angle of rotation was measured at the middle of span by ‘means 
_ Of a reflected light-beam. The speciman used for this test was an aluminum ~ ; 


, alloy beam seventy inches long and having the balanced Z-section shown el : 
‘Fig. 5. ‘Free warping conditions at the supports were attained by using 
e _roller-type bearings oriented so as to prevent both rotation and lateral dis- - 
rs ' _ placement of the supported cross- -section in its plane but allowing complete 
freedom for longitudinal strain. For a load Pi = = 10 lb. applied with eccen- se 
tricity e = 2.603 in. in span 1, the measured value of ¢ in span 2 was -0. 00771 
radians as compared with a calculated value of -0. aehieed radians. 


was applied to a small at the middle of span 1 and the 
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The calculated value of ¢ is founc found ime 1 the following ae 
Aluminum alloy, E = = 10. 5(10)8 psi, psi, G = =4. .0(10) psi. ee. 
40 Length m of center- line of cross- -section (see Fig. 5): 
= 1.319 + + 0.2595 x a+ 0.919 x2= = 3. in. 
= 676)(0. .081)° x 4 = 2605 lb.-in?, 
[ 4a (d+ 4b+2h) 12nd” (b+h) + 2h ?4(2b+3h) + 
= 16284 [ [o. 2.385 + 5.232 + 5. 687] = 224100 Ib.- i 
C and Cy d determined, »k= C/c C,= yo. 01162 = 0. 1078 ir in. and with 
, kf, = ify = 3.773. 


Numerical Cateulations 


+ J0.7357 + 0. 7387] +0=0.3 
from which = 0. 


Elimination of between the first and second of Eqs. (5), after 


hg Sales j- (Pedy Pe ala 


Applied to the a span problem in hand, this becom mes” sin 


ubstitution of the value ra already found above, together with n and vs 


toretio (0. 17352) =. 1(0. -3516) = =- 0. 0000675( Pe) 
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_ Finally, the angle of twist at the mid- — of span 2 is obtained from Eq. 


the , and Ty, found, this gives 

» 


Sham for nee lb. and e = 2.603 in 


3 x 26.03 = -0, 00773 
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SILVERMAN, 1 A. M. ASCE. —The w writer is indebted to Prof. ‘Selberg 
or his comments. W Prof. Selberg questions the validity of Eq. (7). The follow- ¥ 
ing assumptions in the of (7): 
(a) The increase | in the H of inertial forces 
be neglected. This assumption has led to in 

a The bending moment in the suspended span, as the result t of al 


a 


When H is of such magnitude, as for the extreme case of zero sag, so that 


_u may be assumed to be zero, Eq. (7) becomes that of a simply supported span ’ 
ape Hos in a medium furnishing supporting forces proportional to the de- 
flection v. _ Another marge case occurs when the stiffness of the suspended do 
span is zero, i.e. EI = 0. . Eq. _ (7) reduces in that case to that for a simple 


pendulum wt the constamt is BO fat 


Proc. Paper 1292, July, 1957, by I. K. Silverman, 
i. Engr., Bureau of Reclamation, U. S. Dept. of the Interior, Denver, ane 
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= P. P. BIJLAARD, 1 M. ASCE.—Messrs. Haaijer ant Thirlimann, in their 


dealing with buckling of structural steel plates in the horizontal part of the 
cor strain diagram » Say on page 1581-12 that, from the mathematical theo- 
7 ry of plasticity, the stress-strain relations should be of the incremental type | 
and that only the stress-strain law used by Handelman and Prager satisfies 
_ this condition, so that all other theories (including that of the discusser) can + ’ 
be discarded. They refer to tests on combined compression and torsion on * 
_ tubes of mild steel in the horizontal part of the stress-strain diagram at 
~ confirm the stress-strain laws of simple flow theory. They then proceed to 
derive formulas which apply to this special case, for which they mainly refer 


i: ‘It is the purpose of the present discussion to ‘show that the authors’ theo- 
‘retical results, although not stated t that ut way, are actually based on the discus- 
Be | ser’s theory, as published in various papers in the past twenty  years(2, 3,4,5,6, 
7 4, 8) In the first place the discusser’s theory includes both deformation and i. 
4 incremental theory, as pointed out in all his papers and has been stated else- 
: Ww where. (9) The stress-strain relations for simple flow theory as derived by 
Handelman and Prager (1948), follow directly from the discussor’s 
_ published in 1938, by assuming e= = Ep/Ee , = 0, where Ep and Ee are the equiva- 
3 plastic and elastic strains at unelptont buckling and for the special case 
_ of uniaxial compression discussed by them ( = 0) (reference 2,p. 480 and 
Pet, 3, p. 733), or, with other notations, for E, s/E = 1, where Eg is the secant 


modulus at buckling (ref. 6, ‘Vand 8). _In papers published before 


1 ‘much too high theoretical buckling stresses, the discussor assumed defor- 
a. a mation theory, which he showed to give a lower bound for the buckling stress, 


estimate the buckling stress (ref. pp. 133, 734). For »uckling below the 
yield point, the only case that had to be considered: for structural steel, €, 
and hence e is Small » that incremental or f flow theory gives only 


were alreacy much higher than according to Chwalla’s theory ~ 
4 _ the same plastic re reduction factor for plates as for columns) which at that an 
time in Europe was assumed to ‘be right. He also showed that experiments on 
- steel plates by Chase , buckling at the yield point, were in good agreement with 


a a. Proc. Pape Paper 1581 1, April, 1958, by Geerhard Haaijer and Bruno ‘Thirlimann. 
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; mation theory, his theoretical results were in excellent agreement with these 
- experiments and further developed his theory for aircraft structures (ref. 6 . 
: and numerous later papers in the Journal of the Aeronautical Sciences), show- ney 
_ ing also excellent agreement ' with tests by the NACA and others. (8) His theory | 
is now generally used in aircraft design,(11, 12) sometimes with Stowell’ s 
_ simplification (Poisson’s ratio in the elastic range is 0. eo a 
__ That in certain cases deformation theory would lead to impossible results 
: was already shown by the discussor on page 734 of reference 3. However, 
_ this does not mean that in the special case of buckling, where the ratio be a 
>: tween the stress components changes very smoothly, the actual behavior can- “2 
not be similar to what would follow from deformation theory, although it may 
_ be different in other cases than buckling, that is, for non-linear behavior. In 
references 6 and 8 the writer showed how this behavior, which is in ac- © re 
‘cordance with buckling tests on aluminum alloys, can be explained. . It does 
- not conflict with the theory of plasticity. (13) it places the general behavior | 
between that according to deformation and simple flow (incremental) theory, | 
as it actually is. (6,7 1) In this connection it may be noted that, although in this 
- country flow theory is used for problems other than buckling, in Russia defor- 
mation theory is used, Batdorf and Batheneky's s slip theory, which is theo- 
retically sound, but was ; not completely confirmed by non- buckling | testson 
compressed and twisted tubes, for the case of buckling also reduces to defor- _ 
mation theory. (14) ’ The actual behavior of a polycrystalline metal under all an 
a 3 kinds of stress variations could only be ‘predicted by a a theory based on the 
oa : interaction of its dislocations, flaws, etc., which does not exist, but as ex re. 
_ plained in references 6 and 8 B and others, for the special case of buckling me 
ie materials like aluminum, magnesium, titanium alloys and high strength = 
steels, and in an approximate way, for mild steel below the yield point, the - 7 
= behavior is similar to that according to deformation | theory (but may be gener- 
incremental). 
of eccentricities, as suggested by Onat and Drucker, (16) does not lead to an 
agreement between simple flow theory ar and buckling tests. ere 
_ _-However, as already observed on page 55 of reference 17, it is completely — 
> understandable from the discusser’s theory that in the horizontal part of the | 
_ ‘stress- strain diagram of mild steel, which is considered by the authors, the 
_ behavior is close to simple flow (incremental) theory. His explanation of the 5 
behavior of aluminum alloys(6, 8) is paces on the fact that, due to unavoidable ? 
= stress increases ; during buckling. In the horizontal | part of the ~~ 
_ stress- strain diagram this increase is impossible, so that the same case oc- 
curs as in a tube which, after subjecting it to compression, is twisted while a 
keeping the | compressive stress constant. In this case, as follows from the 
es" authors’ tests and also from earlier tests by Hohenemser and Prager, (18) to a 
which the discusser referred in reference 2 and 3 (pp. 480 and 733 , Fespective- 


— ly), the material behaves indeed in accordance with incremental a malt 


_ The authors assume that in this horizontal part of the stress-strain dia- 


gram, at buckling the entire cross section behaves in a plastic way, based on 
a proof given in reference 1. The discusser already has proved that in a oo 
similar way in references 2 and 4, | pages: 416 and 52, respectively, ae 


Hence, here the discusser’s theory applies, withe =QOorEs/E=1. 
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H On page 2 and in their summary the authors say that *til now it was wasas- 
sumed that at the yield stress a plate has no rigidity against buckling. As 
pointed out in reference 2 the discusser actually developed his ‘ theory on the 


Tie buckling of plates with the purpose of showing that at the yield — 7 
a plate ha has still a considerable r rigidity against buckling (as also noted in aa 


f 6, h f f the h’ t). ‘He 
: 7 re erence 6, to prove his t eory on the de ormations of the earth’s crus e 
‘anal the — point ‘in later papers and in a report to the Column Re- 


19) so that this may be to be known since a long 

‘For the basic formulas in the authors’ paper, for which they refer to refer- 

nce 1, they could more appropriately have referred t to poe” 
Papers. In the discusser’ notation their constants are 


= 


x 


and Poisson’s ratio v 


v in the elastic range, where in 1 the authors’ notations, 


can, for “example, be obtained from reference 6 as functions of E E, Et Pi, Po 
Ox and | P2 = example, their 


) an at 
Hence, from Eqs. 6. 4) and (6. 


formute follows directly from Eqs. 


n I (10) and (18) of the discusser’ s paper 
reference 5, which is the authors’ reference 9, which equations also occur 


in an earlier report of the discusser to the sixth International seca of . 
4 Applied Mechanics, Paris, 1946. They give 
where in the authors’ notations ome of 1.8) 


p= = infinite, so > that one ay 


bt 
2). 


In the same way Eqs. 6. 0 and (6.7) of reference — 
1, from which the te Os Dx, Dy and Dy ' are computed, are identical 
the discusser’s formulas(19) in reference 


. As a last Point the discusser may refer to the authors’ Eqs. ( (10) to 


which, instead of deriving them in reference 1, without reference to the dis- Bn 3 
cusser’ ’s results, could have been obtained from Table of reference 
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elastic buckling stresses de have to be multiplied ‘with to obtain plastic 
stresses” Ccr- This would also have avoided errors as occur in 
tl their cases (3) and (4) will be 
(1) Loaded edge es x = 0 and x = Q hing ed, unloaded edge y = 0 hinged and 2 
a ee. _ loaded edge y = b free. . Using Eqs. (1), with v = 0.3, from Table 5 of | 


is identical to Eq. (34) of reference 


2) Loaded edges x = 0 and x = = I hinged, ‘unloaded ¢ edges + d/2 hin id. 
From Table * of reference 6 


ngica 


accordance the authors’ ‘Eq. since Dyx = 
.T (3) Same as case 2, but unloaded e t d/2 fixed. meen 


which is practically in accordance the authors’ (17). 
(4 4) Same as case 1, but unloaded edge 


1, but unloa = 0 fixed. From Table 5 of refer- 
e authors | in (13) 


‘(0.765 (2, D, - +1 40 
whch is is similar to Eq. (3. 15) of reference 1. Itis obvious that t 
authors’ ‘formula (13) must be in error, since in the elastic oot ‘it ae 
does not reduce to the elastic buckling 
Since in the elastic aC, = Dy = E/(1- v2), 
— = v?) ) and ‘ 
v = 0.3, the pire Ed (13) y yields 
Hence one should better use the se the present Eq. (4). 
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the of requires the assumption of un 
known eccentricities for determining the buckling stress or critical strain, ‘it a 
is necessarily very arbitrary. For example, their value 4; = 2400 ksi is 
based on Fig. 10 of reference oh from which Gee ¢ could have been chosen just — 
as well as 1500 or 3000 ksi. . Their method is therefore not capable of predict- 
ing the critical strains but has to be based on tests. The authors should a 
therefore be complimented on their careful and important experiments, on 
which their selection of the values for G t» Dx, Dy and Dxy i is A sanity based. " 
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INCREMENTAL COMPRESSION TEST FOR CEMENT RESEARCH? 


. HRENNIKOFF, 1 ASCE. .—The writer is somewhat disappointed in the 
‘small extent of the response to his paper describing what he believes to be a ~_ 
new approach to the study of creep in cement. This makes him particularly — 
grateful to the sole contributor to the Mr. ‘Keith for his 


aaa s work so far involved no concrete but only cement and partly mortar, 3 
xs feels that the creep of concrete originates primarily in cement, pam i 


incremental tests, namely the regularity of values of the coefficient of lateral 
deformation The e values 0 of during all 
is ; changed are virtually constant and to 0. 2, both in the moist 
cylinders and in the dry ones of the same composition. On the other hand, 
_ the values of 4 observed in the moist cylinders while the load is kept consta ant — 
are different at different load levels, although they tend to remain constant at a 
_ each particular level. While on the way up, the values of yw are close to zero 
on the lower steps of 10 and 15 kips, and they grow from step to ) step rising 
7 high as 0.2 or even higher on the 35 kips step. On the way down the values | 
_ of wu during creep are close to zero on the upper steps and increase to 0.10 - | 
i 0.12 on the > levels 0 and 5 kips. This behaviour strongly suggests that creep _ 
is not the same kind of deformation, only developed gradually, as the one that | 
is taking place while the load is being changed, but is something entirely) 
j “3 The r regularity in the values of the uw stated here has been observed con- __ 
__ sistently in a great number of tests, but the exact magnitude of the lower a 


involved and the difficulty of obtaining the necessary ry sensitivity of the ‘measur- 

ing apparatus. Thus the longitudinal deformation at 10 kips level in Fig. 9 is iz 
estimated at 23 micro-inches per unit. An error of little as one micro-inch _ e 

per inch in reading the lateral deformation would change u by a substantial 2 

of 1/23 - 0.043. Errors of two or three times this amount may be 


‘metal for its segments and isoelastic springs for the sensitive elements. 

24 While during the load changes and on the constant load steps with great amount pi 

Bros. Pape Proc. Paper 1604, April, 1958, by A. Hrennikoff. 
Dept. of Civ. Eng., Univ. of British Columbia, Vancouver , B. c., 
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of creep the inaccuracies | in n the > computed mechanical properties, wm 
_ they are such as to make the computed values of H unreliable. — tae 
_ Most of the tests including the ones involving the specimen C.58 58 (Fig. 9 and 
A 10) show negative values of u on the lowest horizontal step, on the way up, oe 
_ indicating a decrease rather than increase of the diameter of wed specimen, in — 
the course of contraction of its length u under the action ¢ of creep. _ The writer a 
attributed the negative reading of the lateral deformation to the errors of the 
- gauge and assumed the corresponding true value of uw to be zero. There is 4 
however | a possibility that the contraction of the d diameter is a genuine phe- — q 
omenon brought about in the following manner. As” the load on the specimen © 
‘is quickly increased its Length is elastically decreased, and so is the volume, 
as the elastic value of u = 0.2 is less than the value 0.5 ) required fo for the con- < 
— stancy of the volume. The decrease in volume is accompanied by adecrease 
in void space, occupied largely by water, which is thus brought under com- — 
pression, producing a bursting action on the cylinder. The cylinder is thus a) 4 
_ extended laterally not only by the elastic Poisson’s ratio effect, but also. by | 
the tension created by the pressure of water occupying its voids. In the period 
of time when the load is kept constant the water, in an action similar to tae, * 


of the longitudinal ‘stress in n the c cement skeleton, ‘causing its longitudinal 7 
creep, and on the other hand it causes a contraction of cement in the lateral | 
direction, occasioned by the disappearance from it of the lateral tension. bent 
a The effect of the bursting action of water desc ribed here on the longitudinal — 


creep and on the lateral negative creep may be true, but its extent, as is be- __ 
"lieved by | the writer, is very small. ma Instruments of extreme sensitivity a and - 
accuracy, greatly exceeding the ones used by the writer, are needed to prove — 
the existence of negative lateral 
. In his more recent studies the writer has discovered what appears to a 

- the major cause ¢ ofc creep in moist ‘cement. ‘It is the ‘slow - deformation of the — 
water film enveloping the elements of cement gel. The film is held to cement _ 
bya strong attraction, , which transforms the water of the film into a viscous is 
E solid. The writer er does not wish at this — to go ) more extensively into el ay? 


ru "The writer regrets be to answer with certainty the question 


‘The values of longitudinal creep in in these cylinders | are very small, and of the a 
creep—even smaller. Thus the total creep on all seven 


creep at the higher levels compared to the lower ones, no ser Foarice 4 
_ with regard to ) creep in dry cylinders have been observed. D3 x | 


writer expresses hope that in spite of the lack of his 
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An error of one or two units in the lateral strain would result 
An error of one ears that creep in the dry cylinder 
_ by minor internal breakages, 


| 


DISCUSSION 


ers in the field of cement and concrete. ay 
CORRECTIONS. —Page 1604-12, 6th line the rd paragraph should be 
_ Same page, 9th line > of third paragraph should be “on loading”, , rather than A 


pape Page 1724-4, corrections to the pape paper in July 1088 issue of the coun, 
10” is left out from the figure. 
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THE ELASTIC STABILITY OF OF THIN SPHERICAL SHELLS* 


Discussion by Wen Liang Chen 

LIANG CHEN, i J. M. ASCE.—In the determination of critical 
pressure from the expression of total potential energy it was assumed in this 
apes that the parameter C is constant for a given shell in order to reduce 
the computational work. As pointed out by the author in his Sc.D. thesis (1) 
on which this paper is based, experimental results of Kaplan and Fung(2 a, 

_ showed that for a given shell the mode of deflection changes during the loading 
“Process such that in actuality the parameter C is not constant for agiven ~~ 

Shell. The author reasoned that, since his investigation is concerned with ies 
determination of buckling pressures, rather than deflection shapes, it is , 
possible that an incorrect deflection equation sais yield a fairly good - 

_ During the course of research for the writer’s Sc. D. thesis(3) { the critical | 
pressure was determined from the expression of total potential energy used _—° 

by von Willich in such a way that the assumption of constant C was avoided. ; oD :3 
The total potential energy was minimized with respect to both parameters Wo 7 : 

and C. Computation was carried out on the LB.M. 704 computer in the M.LT. 
Computation Center. Fig. 1 shows the comparison 1 of the theoretical results. 
Curves 1 and 2 are results obtained by the method of minimizing the total po- 


_ tential energy with and without the assumption of constant C, respectively. 


_ The agreement between Curves 1 and 2 is fairly good for lambda equal to or 
- less than 5 with a maximum discrepancy of 15%. However, the discrepancy -9 


increases rapidly to 41% as lambda increases from 5 to 6. Thus, the prox- 
-_imity of Curve 1 to the experimental points. in the range 5 < X<8isacci- 
. dental, since the elimination of one assumption in the computation gives re- 
a Also shown in . Fig. 1 are the results of Archer, (4) given by Curve 3 which 
are obtained from perturbation solution of the governing non-linear differ - 
ential equations, and the results of | Weinitschke,\ (5) given by Curve 4 which are 
4 obtained from 5 power ‘series solution of the governing non- -linear ‘differential Sg 
equations. . The dotted portion of Curve 4 is a lower bound to the critical 
pressure given in Reference 5. It can be seen from Fig. 1 that for A< 5. 55 
there | is good agreement between the results obtained from the minimal _ “ ca 
"principle by the writer and from either the power series solution or the per- a 
turbation solution of the non-linear differential equations. For A> 5.5 Curve 
2 rises very - rapidly and is inaccurate probably because the assumed de- “a a 
"flection with respect to which the total potential energy is minimized is not 


a. Proc. Paper 1897 » January, 1959, by Gideon P. R. von Willich. a >) Se 
1. "Graduate Student, Massachusetts Inst. of Mass. 
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An investigation is reported in Reference 3 in which the effect of initial 
geometrical imperfection of the middle surface of shell on the buckling ~~ 
‘pressure is to be quite: significant and could ; account for the scatter of 
’ 


1. Gideon P. von ‘Willich. The Elastic Stability of Thin Spherical ols 
“’ Sc.D. ‘Thesis, Department of Civil and Sanitary Engineering, Massachusetts 


2 A tol nese. od Aa> ti Gt Devi. 


pee 2 © Kaplan and Y. C. Fung. _ A Non- linear Theory of Bending and Buckling © 
of Thin Elastic Shallow Spherical Shells. N.A.C.A. Technical Note 


Wen Liang Chen. Effect of Geometrical on the Elastic 


Buckling of Thin Shallow Spherical Shells. Sc.D. Thesis, Department of 
and Sanitary Engineering, Massachusetts Ins Institute of Technology, 
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Under Uniform Pressure. Quarterly of Applied Mathematics, Vol. XV, 


LXV, 


x J. On the Non-linear Theory of Shallow Spherical Shells. 
iv Journal of the Society for Industrial and Applied Mathematics, Vol. 6, ; 
No. 3, 300-238 
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Discussion by Turgut Sarpkaya 
TURGUT SARPKAYA,’ 1 A. M. ASCE.—Since its inception by Vaschy and re 
later its emphatic presentation and popularization by Buckingham, 
al analysis has enjoyed ; a particular esteem among engineers as well ee 


mathematicians. There ¢ are, to be sure, rather comprehensive treatments of 

the subject some of which are concerned mainly with its philosophical aspects, 
4 some with its mathematical foundations, and some with its practical appli- 
—s cations, However, any iny effort made toward a better determination of the Limi- 
tations of dimensional analysis and toward the simplification of the formation — 

of Pi-terms, such as made by the a authors, is a welcome contribution. on a 4 

The a question has been raised as to the effectivences’ of dimensional analy- 

sis. This writer believes that as a tool alone it is impressive rather than a 

powerful. It should be made quite clear that no useful result could be deduced an 

from dimensional be ager unless it is we by expe experience, , experi- 


on personal eameadaald it is no better and no worse than the researcher’s | 
4 understanding of the particular phenomenon. It should also be made quite <a 
clear that € even the entire analysis is performed meticulously, the research- 


or regione paradox and some other hydrodynamical para-_ 
Sade (1) ) Furthermore , although « one is closer to his goal after re ll 
Considerations, the final and best grouping of variables and the extent of the — 
attained still remains to be determined. Indeed, one is 
 -~. sure until one ha has all the necessary rts amar results whether 
certain Pi-term or terms are worth forming. __ 
In a practical application, once the independent are ¢ determined 
~and the repeating variables are judiciously selected, the application of the Pi- -. 
theorem itself is a a purely mechanical process, no matter which one of the _ 
following methods is used: conventional determinant method, the method of _ 
spectal arrangement of the the variables in the determinant form, and | the method 
selection of repeating ng variables to represent the kinematic, ‘and Jab 


a fourth method, devised by this writer(2) about s six years ago for the for- 7 


Sorte syatent, (four dimensional for heat transfer, electric potential or | 
charge problems). Only | the basic ideas will be presented |here for practical 
demands mathematical rigour may found in the 


Prof. , Eng. Mechanics Dept., Univ. of Nebraska, , Lincoln, Nebr. 
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An axis L » and an ordinate T is taken as shown in Fig. 1. The scale eon 


Each | point (n.L), or (n. ?, ‘distance from the origin represents Ln, and Tn, 
ff 
respectively. n is an ordinary real number and L, and T, are unit quantities. 
Then the following vectorial are self f explanatory: 
AC /AB 
The nature of the axes is such that to divide or to multiply one vector by 
another, an equivalent and parallel vector must be drawn. If a division is to 
be performed, the vector is drawn in the opposite direction; if multiplication, 
in the same direction. It is significant to note that in the new vector system — 4 
_ the ratio of the lengths of parallel vectors is equal to the ratio of their ex- 


ponents. ‘Iti is extremely important also” ton note that each closed d polygon i in 


clockwise. Referring again to Fig. 1, it seen 


One can now examine the system more generally where the operations are © 
made in a three | . dimensional space in which each point - is represented by (L®, 
and each vector by (Lo. M" As will have become evident 
7 * from the foregoing explanations, the properties of the the two dimensional coordi a 
nate system hold true for the three dimensional case as well. olanamib mon 
Having established the fundamental comoapte of the coordinate system one = 
can employ it in the formation of Pi-terms. . If the. meanings length, mass, and 
are assigned to M, and respectively, then the vectors in Fig. 
represent the dimensions of various physical quantities as OA, 


4 
| AB.BO.CA 
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length; OB, area; OC, volume; OD, OE, Kinematic viscosity, ¢ circe- 
lation, or unit discharge; OF, 
per unit time; OK, surface tension; ON, pressure intensity, elastic modulus, 
‘or shear intensity; OP, specific weight; and OR, the mass density. - Obviously, 
geometric quantities are on L®, axis, kinematic quantities are in(L®, T"), 
plane, and dynamic quantities are in (LA, T", M®), space. Each vector in this 
system represents the dimensions of products which may or may not have a 
physical meaning. . Every: vector such as OF in Fig. 2, actually represents the 
_ n-th power of the dimension vector, (OT = V2), and usually in the case of _ 
n= 1, has a physical meaning. As far as the selection of repeating variables 


: 7 is concerned, it is well known that the determinant of the exponents of repeat-_ 
s variables must not be equal to zero. It can be seen very easily that in the _ 
is 


-vectorial system this mathematical c condition on corresponds to the condition 
1 = that the vectors representing the dimensions of repeating variables ‘must not 
om form a closed polygon. For example, the vectors corresponding to area, ve- - 
locity, and 1 mass density c cannot form a closed polygon r of value 
oftheir exponents, 

- To illustrate the method, the following examples a are given: ee: 


1. The combination of V (v elocity), D (length), and g os 
ation), is desired. In Fig. 3, these vectors are first drawn, and then .. 
_ Closed polygon is formed. — From the small triangle, rotating in the o 
clockwise direction, one obtains: ‘ V.D-1/2.g-1/2, Similarly, from the 


‘The combination of p v (velocity), , and p (pressure), is shown 


> 
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(force), A density. The numbers 


The combination of u (dynamic viscosity), density, velocity, and length 


is shown in Fig. 5; and of surface tension, density, — and length» 


_ The examples given above clearly show that all the other classical di- 
mensionless numbers can be obtained very easily. All is necessary is a little | aa 
‘practice. The application of the method to the systems which have more than ~ 
three | basic units, ie the same. To prove the s a 
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DISCUSSION 
og 


A 


_— involving | the transference of heat is given. The temperature is regarded as ik 
a quantity not yet reducible in terms of the fundamental units of mass, length, © 
_ and time. In Fig. 7, the combination of € (specific heat per unit mass), a 
(heat transferred per unit time), @ (absolute (a 
In the examples given above, axes were taken as MM, and 
‘perfectly permissible t to use Ln, Fa, TA, and 6", which may be more “Coane 
’ 3 to the « engineer who is ‘accustomed to using force rather than mass. - In 


0A. 


ai 


han 


b 


_ electrical engineering problems the axes could be chosen either as L®, Ma, 
§ T, and Qn, or as L®, M®, T®, and $2, in which Q, and ¢, represent 
electric charge and the electric potential respectively, 
in conclusion, it it is hoped tha that the addition of a vectorial method for the a a 
r toumnatign of Pi- terms to those proposed by the authors will add further clari cit 
1s ty and physical understanding to dimensional analysis, and consequently to the 
_ determination of functional equations in the diversified fields of science. = 
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